Quantum computation with unknown parameters 
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We show how it is possible to reahze quantum computations on a system in which most of the 
parameters are practically unknown. We illustrate our results with a novel implementation of a 
quantum computer by means of bosonic atoms in an optical lattice. In particular we show how a 
universal set of gates can be carried out even if the number of atoms per site is uncertain. 
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Scalable quantum computation requires the implemen- 
tation of quantum gates with a very high fidelity. This 
implies that the parameters describing the physical sys- 
tem on which the gates are performed have to be con- 
trolled with a very high precision, something which it 
is very hard to achieve in practice. In fact, in sev- 
eral systems only very few parameters can be very well 
controlled, whereas other posses larger uncertainties. 
These uncertainties may prevent current experiments 
from reaching the threshold of fault tolerant quantum 
computation (i.e. gate fidelities of the order of 0.9999 
lH), that is, the possibility of building a scalable quantum 
computer. For example, in quantum computers based on 
trapped ions or neutral atoms the relative phase of 
the lasers driving a Raman transition can be controlled 
very precisely, whereas the corresponding Rabi frequency 

has a larger uncertainty Afl. If we denote by T the 
time required to execute a gate (of the order of 
then a high gate fidelity requires TAil <C 1 (equivalently, 
Afl/n ^ 1), which may very hard to achieve, at least to 
reach the above mentioned threshold. 

In this letter we show how to achieve a very high gate 
fidelity even when most of the parameters describing the 
system cannot be adjusted to precise values. Our method 
is based on the technique of adiabatic passage, combined 
with some of the ideas of quantum control theory. We 
will illustrate our method with a novel implementation 
of quantum computing using atoms confined in optical 
lattices 0. If the number of atoms in each of the po- 
tential well is uncertain, which is one of the problems 
with this kind of experiments, most of the parameters 
will have an uncertainty of the order of the parameter it- 
self (e.g. |Ar2| ~ Q), which under normal circumstances 
will give rise to very poor fidelities and even impede the 
performance of quantum gates. As we will show, using 
our method not only quantum computation is possible 
but even very high fidelities could be achieved. 

The use of adiabatic passage techniques to implement 
quantum gates is, of course, not a new idea. In fact, 
several methods to perform certain quantum gates using 
Berry phases have been put forward recently |p|, p|, p^. 
Furthermore, adiabatic passage techniques have been 
proposed as a way of implementing a universal set of 
holonomies [ p^ , i.e. quantum gates which are carried 
out by varying certain parameters and whose outcome 



only depends on geometrical properties of the paths in 
parameter space In all these proposals, physical im- 
plementations of standard quantum computation have 
been adapted so that the quantum gates are performed 
in an adiabatic way giving rise to holonomies. Despite 
its clear fundamental interest, it is not clear yet if such a 
novel way of implementing the quantum gates may offer 
real benefits with respect to the original proposals. In 
contrast, in our illustrative example adiabatic passage is 
required to perform quantum gates and therefore it is an 
essential tool not only to achieve the desired precision 
but also to build a quantum computer at all. 

The outline of the paper is as follows. First we will 
show how to produce a universal set of gates (Hadamard, 
phase, and CNOT) starting from two Hamiltonians in 
which only one parameter is precisely controlled. Then 
we will show that this method eliminates an important 
obstacle in a particular physical scenario that has been 
proposed for quantum computation [Q, namely a set of 
atoms in optical lattices interacting via cold collisions. 

Let us consider a set of qubits that can be manipulated 
according to the single qubit Hamiltonian 



and the two-qubit Hamiltonian 

i?2 = A|11)(11| + ^I® (a+e''^ + 
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As mentioned above, we will assume that most of the pa- 
rameters appearing in these Hamiltonians are basically 
unknown. On the other hand, we will not consider any 
randomness in these parameters (the corresponding er- 
rors could be corrected with standard error correction 
methods as long as they are small). In particular we 
will assume that only ip can be precisely controlled. For 
the other parameters we will assume that: (1) they are 
given by an unknown (single valued) function of some 
experimentally controllable parameters; (2) they can be 
set to zero. For example, we may have n — /(/), where 
/ is a parameter that can be experimentally controlled, 
and we only know about / that /(O) =0 and that we can 
reach some value flm = film) 7^ for some Im 0- Below 
we analyze a particular physical scenario which exactly 
corresponds to this situation. However, we want to stress 
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FIG. 1: Schema of how the parameters of Hamihonian (J^ 
have to be changed in order to perform a phase gate (a), and 
Hadamard gate (b). In Fig. (c) we show a possible evolution 
of the parameters A and O for the Hadamard gate. 



that this situation can be naturally found in more gen- 
eral scenarios. For example, for atomic qubits the states 
|0) and |1) may correspond to two degenerate atomic 
(ground state) levels which are driven by two lasers of 
the same frequency and different polarization. The cor- 
responding Hamiltonian is given by (|l|), where the pa- 
rameters (fi, f2, A describe the relative phase of the lasers, 
the Rabi frequency and detuning of the two-photon Ra- 
man transition, respectively. The Rabi frequency can be 
changed by adjusting the intensity of the lasers, and the 
detuning and the phase difference by using appropriate 
modulators. In practice, f2 (A) can be set to zero very 
precisely by switching off the lasers (modulators) and (p 
can be very precisely controlled to any number between 
and 2tt. However, fixing 51 or A to a precise value (for 
example, 23.098 kHz) may be much more challenging. 

The idea of obtaining perfect gates with unknown pa- 
rameters relies on the combination of adiabatic passage 
techniques [|| and ideas of quantum control Let us 
recall the basic idea of adiabatic passage. Suppose we 
have a Hamiltonian that depends parametrically on a set 
of parameters, denoted by p, which are changed adiabati- 
cally with time along a given trajectory p(i). After a time 
T, the unitary operator corresponding to the evolution is 
given by 



$„[p(T)])(ci>„[p(0)]|. (3) 



Here, \^a{p)) are the eigenstates of the Hamiltonian for 
which the parameters take on the values p. The phase 
(pa is a dynamical phase that explicitly depends on how 
the parameters p are changed with time, whereas the 
phase ipa is a purely geometrical phase ands depends on 
the trajectory described in the parameter space. Our 
basic idea to perform any given gate is first to design 
the change of the parameters in the Hamiltonians (|l|)-(0) 



such that the eigenvectors evolve according to the desired 
gate, and then to repeat the procedure changing the pa- 
rameters appropriately in order to cancel the geometric 
and dynamical phases. 

Let us first show how to perform the phase gate U = 
^10(7^/2^ We set A for all times. The parameters 
(ri. If) have to be changed as follows [see Fig. ^a)] : 
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Steps (i,ii) and (iv,v) are performed adiabatically and 
require a time T. The double arrow of step (iii) indi- 
cates a sudden change of parameters. Note that 57(0) = 
n{2T) = 0, n{t) = n{2T-t) and (p{t) = Tr + 0-ip{2T-t), 
which does not require the knowledge of the function / 
but implies a precise control of the phase. A simple anal- 
ysis shows that (i-v) achieve the desired transformation 
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1). Note also that the dy- 



namical and geometrical phases acquired in the adiabatic 
processes (i-v) cancel out. 

The Hadamard gate can be performed in a similar fash- 
ion. In the space of [A, il^ = cos{ip)], the protocol is 

(0,r!„0 ^ (A„,a„0*^^A„„0)*^^A,„,5]™)'^^ 

(0, a„) ^ (0, -r!„0 (A, -a„) '^'^ (A, 0), (5) 

as shown in Fig. ^(b-c). In order to avoid the dynamical 
phases, we have to make sure that steps (i-v) are run in 
half the time as (vi-vii). More precisely, ii t < T, we 
must ensure that A{t) ^ A{T - t), fl^{t) = n^{T -t), 
A(T + t) = A{t/2) and n^{T + t) = n^{t/2). With this 
requisite we get ^(|0) + |1)) ^ |0), ^(|0) - |1)) ^ 
— Again, the whole procedure does not require us to 
know 57 or A, but rather to control the evolution of the 
experimental parameters which determine them. 

The C-NOT gate requires the combination of two 
two-qubit processes using H2 and one local gate. The 
first process involves changing the parameters [A, f2a: = 
Cl cos(iy9)] of Eq. (^ according to 

(A„„0) ^ (A„,57„) ^-^^ (0,an) ^S^^0,-f7„) 

(A„,-an) ^ (A„„0). (6) 
This procedure gives rise to the transformation 

Ui = \0){0\(g)l+e'^\l){l\(g)icry, (7) 

where ^ — 6{t)dt is an unknown dynamical phase. 
The second operation required is a NOT on the first qubit 
U2 = (|0)(1| + |1)(0|) m. Finally, if A(i)(t) denotes the 
evolution of A in Eq. we need to follow a path such 
that A(3)(t) = A(i)(i), n^^^{t) = 0. If the timing is 
correct, we achieve the transformation 



C/3 = (|0)(0|+e'?|l)(l|)®£ 



(8) 
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FIG. 2: Log-log plot of the gate error, E = 1 — T, for 
different local and nonlocal gates, (a) Realization based on 
Hamiltonians (|l|)-(|^, with fixed parameters and varying time. 
We plot the error for the Hadamard (solid), phase (dashed) 
and CNOT gates (dotted), (b-d) Errors for our proposal with 
atomic ensembles. We fix the the total time all parameters 
except the interaction constant Ubb- In (b-c) we plot realiza- 
tions for n = 1 (solid) and n = 3, 5 (dashed) atoms per lattice 
site. In (d) we plot simulations with a difference of atoms of 
|n — m| = 0, 1, 2 between both wells (solid, dash and dotted 
lines). The vaues of the parameters are given in the text. 



Everything combined gives us the CNOT up to a global 
unimportant phase Ucnot = |0)(0| ®I+ ® ic^y = 

e-'«C/2C/3C/2C/i. 

In Fig. ||(a) we illustrate the performance of our 
method, as well as its sensitivity against non-adiabatic 
processes. As a figure of merit we have chosen the gate 
fidelity \Tr{ul^^iUreai}?ld'', where d is the di- 

mensionality of the space (2 for local gates, 4 for two- 
qubit gates), U ideal is the gate that we wish to pro- 
duce and Ureal is the actual operation performed. As 
expected, for fixed parameters {Am/rim = ^m/^m — 
1/10, ipm — ■'i'/4}, the adiabatic theorem applies when the 
processes are performed with a sufficiently slow speed. 
Typically a time T ~ 300/rim, 300/r2m is required for 
the desired fidelity T ^ 1 — 10"''. 

Let us now consider a set of bosonic atoms confined in 
a periodic optical lattice at sufficiently low temperature 
(such that only the first Bloch band is occupied). The 
atoms have two relevant internal (ground) levels, \a) and 
1 6), in which the qubit is stored. This set-up has been 
considered in Ref. where it has been shown how single 
quantum gates can be realized using lasers and two-qubit 
gates by displacing the atoms that are in a particular in- 
ternal state to the next neighbor location. The basic 
ingredients of such a proposal have been recently real- 
ized experimentally [p2|. However, in this and all other 



schemes so far Q it is assumed that there is a single 
atom per lattice site since otherwise even the concept of 
qubit is no longer valid. In present experiments, in which 
the optical lattice is loaded with a Bose-Einstein conden- 
sate 1^, this is not the case (since zero temperature 
is required and the number of atoms must be identical 
to the number of lattice sites). We show now a novel 
implementation in which, with the help of the methods 
presented above, one overcomes this problem. 

For us a qubit will be formed by an aggregate of atoms 
at some lattice site. The number of atoms forming each 
qubit is completely unknown. The only requirement is 
that there is at least one atom per site 0. We will 
denote by rik the number of atoms in the fc-th well and 
identify the states of the corresponding qubit as 



\0)k = 



,t("fc-l) 



\l)k=bl 



(9) 



where Oj, (5jl) are the creation operators for one atom in 
levels I a) and |6), respectively. The quantum gates will 
be realized using lasers, switching the tunneling between 
neighboring sites, and using the atom-atom interaction. 

In the absence of any external field, the Hamiltonian 
describing our system is 

H = - jr Ei^Uh + bib,^,) +^y: f^'A^'^b'^- 

k k 

(10) 

Here, Ubb and J^^^ describe the interactions between and 
the tunneling of atoms in state \b). We will assume that 
J^*"' can be set to zero and increased by adjusting the in- 
tensities of the lasers which create the optical lattice. We 
have assumed that the atoms in state \a) do not interact 
at all and do not hop, something which may be achieved 
by tuning the scattering lengths and the optical lattice. 
Both restrictions will be relaxed later on. The Hamilto- 
nian (^h possesses a very important property when all 



J, 



(&) 



0, namely it has no effect on the computational 



basis (i.e. H\'9) = for all states |^) in the Hilbert space 
generated by the qubits). Otherwise, it would produce a 
non-trivial evolution that would spoil the computation. 

We show now how a single qubit gate on qubit k can be 
realized using lasers. First, during the whole operation 



we set J, 



(b) 



in order to avoid hopping. The laser 



interaction is described by the Hamiltonian 



(11) 



For Ubb ^ |Afc|,|r2fc|, we can project the total Hamil- 
tonian acting on site k by an effective Hamiltonian act- 
ing on the qubit which resembles (|l|), with A = A^ and 
f2 = VLk^/uk — 1. Thus, using the methods exposed above 
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we can achieve the Hadamard and phase gates with a high 
precision, even though the coupUng between the bosonic 
ensemble and hght depends on the number of atoms. 

For the reahzation of the two-qubit Hamiltonian (|^) 
we need to combine several elements. First of all we need 
the Raman coupling of Eq. (|l^) to operate on the second 
well. Second, we need to tilt the lattice using a magnetic 
field [Q which couples to states \a) and \b) differently, 
Huh = J2k kgi'^k'^k — bkbk)- And finally we must allow 
virtual hopping of atoms of the type \b) {J^''^ <C \Ubb — 
g\). After adiabatic elimination we find that the effective 
Hamiltonian depends on the number of particles in the 
second site, n2. 



H. 



eff _ 



g - Ubb 



\ll){ll\ + y/n2-ini(g)a^. (12) 



The identification with Eq. (|^) is evident, and once more 
the use of adiabatic passage will produce gates which are 
independent of the number of particles. 

We have studied the different sources of error which 
may affect our proposal: (i) Ubb is finite and (ii) atoms in 
state I a) may hop and interact. These last phenomena are 
described by additional contributions to Eq. (pOl) which 



are of the form Jl°'\al.ak+l+a].^■^^ak), Uaaa\a\akak, and 
Uabo\.bl,akbk- The consequences of both imperfections 
are: (i) more than one atom per well can be excited, (ii) 
the occupation numbers may change due to hopping of 
atoms and (iii) by means of virtual transitions the effec- 
tive Hamiltonian differs from (^ and (||). The effects (i)- 

(ii) are eliminated if (n/Ubb)^ < 1 and {j]^^ /Ubbf < 1- 
Once these conditions are met, we may analyze the re- 
maining errors with a perturbative study of the Hamil- 
tonians ( pl| ) and ( p^ plus the terms (J''^\ J7ab, Uaa) that 
we did not consider before. In Eq. (|l^), the virtual ex- 
citation of two atoms increments the parameter A by an 
unknown amount, Ae// A + 2r2^n/c/(A -|- Uab — Ubb)- 
If Uab ^ Ubb and VL^UkT /Ubb ^ 1, this shift may be ne- 
glected. In the two-qubit gates the energy shifts are in- 
stead due to virtual hopping of all types of atoms. They 
are of the order of max(j(^\ J^°'''Y / g"^ ~ J'^ /Ubb, and for 
J'^T/Ubb ^ 1 they also may be neglected. 

To quantitatively determine the infiuence of these er- 
rors we have simulated the evolution of two atomic en- 
sembles with an effective Hamiltonian which results of 
applying second order perturbation theory to Eq. (10), 



and which takes into account all important processes. 
The results are shown in Figs. ||(c-d). For the two- 
qubit gate we have assumed Uaa = Uab, J^"''' = J^'''' = J, 
Jm = O.OSC/ab, Vtm = J^/IO, g ^ Ubb + Uab/2 and opera- 
tion time, T — 200/r2,„, while changing the ratio Ubb/Uab 
and the populations of the wells. For the local gates we 
have assumed A^, — 1, Q.ra — C/&f,/10 and different occu- 
pation numbers Uk, and we have also changed Ubb/^m- 

We extract several conclusions. First, the stronger the 
interaction between atoms in state |6), the smaller the 



energy shifts. Typically, a ratio Ubb = 10'*f/ab is required 
to make J- = 1 — 10~^. Second, the larger the number 
of atoms per well, the poorer the fidelity of the local 
gates [Figs. ||(b-c)]. And finally, the fidelity of the two- 
qubit gate presents a small dependence on the population 
imbalance between wells. 

In this work we have shown that it is possible to per- 
form quantum computation even when the constants in 
the governing Hamiltonians are unknown. We have de- 
veloped a scheme based on performing adiabatic passage 
with one-qubit (|l|) and two-qubit (||) Hamiltonians. With 
selected paths and appropriate timing, it is possible to 
perform a universal set of gates (Hadamard gate, phase 
gate and a CNOT). These procedures cannot only be 
used for quantum computing but also for quantum sim- 
ulation jTsj. Finally, based on the preceding ideas, we 
have proposed a scheme for quantum computing with 
cold atoms in a tunable optical lattice. Our scheme works 
even when the number of atoms per lattice site is uncer- 
tain. Note that, these ideas also apply to some other 
setups like the microtraps recently realized in Rcf . |l6| . 

Wc thank D. Licbfried and P. ZoUcr for discussions and 
the EU project EQUIP (contract IST-1999-11053). 
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